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Abstract. Jeans instability of finite massive bodies at hydrostatic equilibrium is stud- 
ied. Differential equation governing the evolution of infinitesimal disturbances is derived. 
We take into account radial inhomogeneity of mass density and other fluid parameters at 
the equilibrium state. Dispersion relation and a simple analytical formula, generalizing 



' the Jeans criterion of instability, are derived 

o ' 



1. INTRODUCTION 

It is well known (See, for instance pQ - that the Jeans' instability forms the basis 

of our understanding of gravitational condensation. In particular, Jeans' mass criterion 
> 

is invoked in astrophysical theories of the formation of stars, gaseous clouds, etc. Usually 

^ ! gravitational instability is analysed in terms of the Jeans' wavelength jSj 
O. 

S: A J =S (1.1) 



Gp 



O 

or, equivalently, in terms of Jeans' mass Mj ~ Po-^j- in this formula G is the gravitational 
constant, p is the unperturbed mass density and c s is the adiabatic sound speed. As is 
now widely known, perturbations in homogeneous fluid with mass greater than a critical 
value Mj may grow producing gravitationally bounded structures. In the process of 
^ . their evolution this structures can achieve states of hydrodynamic equilibrium like stars 

polytropes or gas clouds when pressure gradient equals gravitational force. 

In this paper we investigate hydrodynamic equilibrium and stability of finite self- 
gravitating fluid mass with inhomogeneous distribution of mass density, pressure and 
temperature along the radius. In linear approach we get Schrodiger-like |H| equation wich 
eigenvalues and eigenf unctions give us increments and profiles of disturbances. 

2. BASIC FORMALISM 

Consider spherically symmetric fluid body with radius R and mass M. We assume that 
the system is non-rotating and non-expanding. It can be star, gaseous cloud, etc. The 
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evolution of a self-gravitating fluid is described by the conservation equations for mass, 
momentum and specific entropy, coupled with the Poisson equation 

% + Vpv = (2.1) 

p^ + p(W)v = - Vp - pV<p (2.2) 
- + (vV) S = (2.3) 

VV = 4vrGp, (2.4) 

where v is the velocity, p is the pressure, p is the mass density, s is the specific entropy and 
ip is the gravitational potential. The linearization procedure has clearly invoked that local 
state variables deviate from their equilibrium values through linear fluctuation, namely 

P =Po(r) +Pi(t,r) , p = p (r) + pi(t,r) , 

(2.5) 

s =s (r) + si(t, r) , <p = (p (r) + <pi(t, r) . 

Velocity v itself is infinitesimal. Substitution of ()2.5|) in ()2.1|) - ()2.4|) constitutes equations 
for equilibrium state (P) 

^Vp (r)-V<p (r) = (2.6) 

Po{r) 

VVo(r) = 4vrGp (r) (2.7) 

and for perturbed parameters 

+ Vpo v = (2.8) 

Po-t^ = -Vpi - piVtpo - PoV(pi (2.9) 

^ + vVs = (2.10) 

VVi=47rGp 1 . (2.11) 

These equations must be coupled with the proper equation of state. To simplify the 
analysis and exclude buoyancy forces (|2j,[in]) we take it to be adiabatic 

P~P 7 , (2.12) 
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where 7 is the adiabatic exponent. Thus we get no entropy desturbances 

s(t,r) = so(r) = const (2-13) 

and pressure and mass density disturbances are bound together with adiabatic equation 

Pi = c 2 sPl . (2.14) 

Equations (j2.6J) . ()2.7|) and (j2.12|) for equilibrium state are resolved with well known Emden 
functions (pQ) for polytropic model. For further references we summarize here basic results 
of polytropic theory of ideal gas when polytrope exponent n — 1/(7 — 1). In this case 
radial dependancy of equilibrium parameters is 

Po (r)=p(o)er 1 (0 

p (r) = p(0)Q n n (O (2.15) 
c s 2 (r) = c 2 (0)9 n (0 , 

where £ = r/R and n are the non-dimensional radius and temperature respectively. To 
satisfy equilibrium equations (j2.6|) . ()2.7|) and boundary conditions O(0) = 1, 0(1) = 
parameter 

& = (2.16) 

nc 2 s {0) 

must have unique value for each n. For example for n = 3/2 it equals to £ 3 / 2 = 3.65 and 
for n = 5/2 £5/2 = 5.36 (P). 

Equations ([2.8)1 . ()2.9|) . ()2.11|) and ()2. 14|) with proper boundary conditions forms the 
basis of our treatment of stability. 

Boundary conditions we discuss later. But now we transform perturbed equations into 
the more convenient form. Substituting ()2.11|) into (|2.8|) gives us 

'v| + J=0. (2.17) 



AnG dt 

Taking into account that div rot\I/ = we get from (|2.17j) 



^V^ + rot* = 0, (2.18) 



4ttG dt 

where ^ stands for vector-potential of the flow p v. As we can see from (12. 18(1 . vector- 
potential ^ represents the "axial part" of desturbances and is not derectly bound with 
gravitational potential. So we assume it to be zero. Next, inserting ()2.18|) with St = 
and d2H| into fljESJ), we obtain 

1 d 2 

-4^ V -Jr + PoV<Pi = -Vcfo - Pl V<p • (2.19) 
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A harmonic time dependence ~ exp(— iut) of the perturbations can now be assumed in 
terms of constant complex frequency uj so that equation (j2.19j) after simple transforma- 
tions becomes 

V yi = -4,G V ti^' Vy ° . (2.20) 

Taking divergence of equation ()2.20|) and using (|2.11|) we finally get equation for mass 
density disturbance pi only 

Before discussing equation 1)2.21)1 we make the last simplification by introducing auxiliary 
function 

* (r) = I^ (w(r) + c ' (r)) (2 ' 22) 

and replace pi with new function u 

pi = exp(-x)w. (2.23) 
It yields final equation for unknown function u 

exp(x) V ( e ^~ X ]f s \7u) + u = 0. (2.24) 



u 2 + AnGp 

Equation ()2.24j) represents a modification of Sturm-Liouville eigenvalue problem and, 
together with suitable boundary conditions, provides the eigenvalues and eigenfunctions 
for the perturbations. Boundary conditions are summarized as follows. 

On the free moving surface of the body pressure must be equal to zero 

Pl (t,r = R) = 0. (2.25) 

Boundary conditions are stated at unperturbed surface so as small deviations of size leads 
to second order magnitude of perturbance. If the local sound velocity at the body surface 
is not equal zero, Eq. ()2.14|) . ()2.25)) constitutes that mass density is equal to zero 

Pl (t,r = R) = 0, if Cs (r = R)=£0. (2.26) 

But if, as it is in polytropic models, c s (r = R) = then boundary condition ()2.25|) is 
valid with arbitrary p\{r = R). To find out baundary condition for mass density in this 
case we integrate f|2.8|) on the unperturbed volume. Taking into account Gauss theorem 
we get 

j[ Pl dV+ I Po vdS = 0, (2.27) 
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where dS is body surface element. For zero mass density po at r = R second item in 
(|2.27jl is equal to zero and mass conservation low (J2.27)) reduces to 

P!dV = 0, if Cs (r = R) = 0. (2.28) 

A simplification can now be introduced to handle the non-trivial angular dependence 
in the perturbation (J2.24)) . It can be rewritten after decomposing variables into spherical 
harmonics Y/ m 

u{v) = Y J Ui m {r)Y lm {6,<P). (2.29) 

l,m 

Usually perturbations with quantum number I > have small increments so we take into 
account only spherically symmetric perturbations with / = 0. For this case main equation 
in spherical system of coordinates is 

Lu = * U *XP(-X)fs *L) +u = . (2 .30) 

r 2 dr \ u 2 + 4:7iGp dr J 

Operator L has some important features. First of all we now prove that it may have 
only real eigenvalues uj 2 . Let multiply (J2.30)) with complex conjugate u* and function 
r 2 exp(— x) and integrate along radius. We get 

R exp(-x)c 2 du 2 2 f R 2 2> n fo oi\ 

1 r dr + / exp(— x) \u\ r dr = 0. (^-<jl) 



/ uj 2 + AitGpo dr Jo 

First item in 1)2.31)1 is derived with integrating by parts and using that c 2 u equals zero 
on body surface. Subtracting from (j2.31)) its complex conjugate we get after simple 
transformation 

(^ 2 -- 2 *) r , ? P( 7^,J ?| 2 ^ = 0- (2-32) 



jq \lu 2 + AnGp \ 2 dr 

Obviously integral in equation 1)2.32)1 is greater zero and imaginary part of u 2 must be 
zero: Im u 2 = 0. So there are two types of oscillation modes: if u 2 > eigenvalues u are 
real and introduces sound-like oscillations, but if u 2 < eigenvalues are pure imaginary 
and branch with plus imaginary part introduces monotonously growth of perturbation. 

Next, again multiplying equation ()2.3())1 with r 2 exp(— x) and integrating it along r we 
get 

/ 2 zM-x)cl du\ R + r R )ur 2 dr = (2 33) 

V uj 2 + AnGpo dr J Q J 
If c s (R) = first item in ()2.33)) equals zero and equation ()2.30)) automatically conserve 
total mass (cf. I2~2%1 EZHj) 

rR 

exp(— x)wr 2 dr = . (2.34) 



o 



To find out conditions for existance of instability suppose that uj 2 < and mass density 
Po(r) and sound velosity c s (r) monotonously decrease from theirs maxima at r = to 
zero at r = R. Introduce also function of wave number 

k(r) = ^uo 2 + 4TrGp (r)/c s (r) = ^AnG Po (r) - \uj\ 2 /c s (r) . (2.35) 
In this notation equation (|2.30J) can be rewritten as 

exp( X ) d ( r 2*M-x)du\ +u = Q (2 36) 



r 2 dr \ k 2 dr / 
Wave number k(r) equals zero at the radius r where 

AirGpo(r ) - \cu\ 2 = 0. (2.37) 

At this point equation ()2.36|) has singularity and eigenfunction u must have zero derivative 

^(r = r )=0. (2.38) 

Now we build up eigenf unctions of (j2.3(ij) with method like quasiclassic approach in 
quantum mechanics ([8J. "Sewing condition" at the point r = r will give us dispersion 
relation for unstable modes. 

Quasiclassic eigenfunction in the region r < tq finite at r = is 

u = eXp(x/2)fcV2 sin ( ^ k{r) dr] r < r , k 2 > 0. (2.39) 



r \jo / 

As we can see from ()2.39j) . in this region perturbation oscillate with radius. But in the 
region r > r it exponentially decrease 

u = gex P(*/ 2 )l fc l 1/2 exp ^_ jf' 1^(^)1 dr ^j r » r0; p < o, (2.40) 

where C is constant. We will get "sewing condition" from analitical continuation of 
expression ()2.40j) into region r < r through up and down halfplanes of the complex 
variable r — r (|8J). This procedure yields 

Cexpmi&l 1 / 2 / T, /XI \ Cexpd)^ 1 / 2 / f ro /x 3vr\ , 
y )^ J{ 1 ex p(^-y \k(r)\drj — ► ^ sin^y k(r) dr + ^-j . (2.41) 

Right item of (j2.41j) must be equal to right item of (j2.39j) 



Csin( / k(r)dr + ^ J = sin k(r) dr ) . (2.42) 



Writing integral at right part of ()2.42|) in the form 



pro pro 

k{r) dr = k{r) dr — I k{r) dr (2.43) 

1 1 Jo Jr 
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we get equality condition 

7T 

k(r) dr = — h nn 
4 

n = 0,1,2,... ( 2 - 44 ) 
C= (-l) n . 

Expressions (I2.44J1 define discrete set of unstable modes with encrements u n (u 2 < 0) 
and eigenfunctions of type (|2.39p . ()2.40|) . Eigenfunction u n has strictly n zero nodes at 
the interval < r < R. But the mode with n = can not exist as it does not change sign 
and can not obey the mass conservation low (|2.34jl . 

Finally dispersion relation for unstable modes takes the form 

y/4irGpo{r) - \uj n \ 2 vr 

— ^— dr = — \-nn, n — 1,2,3.... 2.45 

c s {r) 4 

Simple criterion of instability immidiently follows from expression (|2.45jl . If we neglect 
\uj n \ 2 and expand integration from r to R then we obviously get for n = 1 condition for 
existance of instability in the form 



o 



»4ff^*». (2.46) 



3. CONCLUSION 



Now we can apply criterion (|2.46J) to models of massive gaseous clouds mainly con- 
sisting of molecular hydrogen. At temperatures < 90K° rotational degrees of freedom 
are degenerated (0) and it behaves as monoatomic ideal gas with adiabatic exponent 
7 = 5/3 and polytrope exponent n = 3/2. At temperatures > 90K° adiabatic exponent 
has standard value 7 = 7/5 and n = 5/2. Criterion (|2.4fij) may be rewritten in terms of 
polytrope model as follows (cf. I2.15[ 12.16)) 



a,. 



Numeric integration in expression ()3.1|) gives for 7 = 5/3 

a 3/2 = 0.92, (3.2) 

and for 7 = 7/5 

Qt 5 / 2 = 1.1. (3.3) 

This result seems to be very intristing. It follows from equations ()3.2)) and ()3.3)) that 
cold hydrogen cloud may be stable relatively gravitational instability, but if the mean 



8 



temperature of hydrogen cloud is high enough, it may be the subject of gravitational 
instability. 

We can also remark here that our result for 7 = 5/3 is applicable to degenerated 
electron gas (white dwarfs) and our proof of its stability is in accordance with well known 
results of the theory of white dwarfs (P,j2],IZj)- 
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